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Abstract  Transmit  beamforming  with  receive  combining  is  a  simple  approach 
to  exploiting  the  significant  diversity  provided  by  multiple-input  multiple-output 
(MIMO)  systems,  and  this  technique  can  be  easily  extended  to  frequency  selective 
MIMO  channels  by  employing  orthogonal  frequency  division  multiplexing  (OFDM). 
Optimal  beamforming  requires  channel  state  information  in  the  form  of  the 
beamforming  vectors  corresponding  to  all  the  OFDM  subcarriers.  When  the  uplink 
and  downlink  channels  are  not  reciprocal,  this  information  must  be  conveyed  back 
to  the  transmitter. 

To  reduce  the  amount  of  feedback  information,  a  new  approach  to  transmit 
beamforming  is  proposed  that  combines  limited  feedback  and  beamformer 
interpolation.  Because  the  length  of  the  OFDM  cyclic  prefix  is  designed  to  be  much 
less  than  the  number  of  subcarriers  to  increase  spectral  efficiency,  the  neighboring 
subchannels  of  a  MIMO-OFDM  system  are  substantially  correlated.  Thus,  the 
beamforming  vectors  determined  by  the  subchannels  are  also  significantly 
correlated.  To  reduce  the  feedback  information  using  the  correlation  between 
beamforming  vectors,  the  receiver  of  the  proposed  scheme  sends  back  only  a 
fraction  of  information  about  the  optimal  beamforming  vectors  to  the  transmitter. 
Then  the  transmitter  evaluates  the  beamforming  vectors  for  all  subcarriers  through 
interpolation  of  the  conveyed  beamforming  vectors.  Since  a  beamforming  vector  is 
phase  invariant  and  has  unit  norm,  a  new  spherical  linear  interpolator  is  proposed 
that  exploits  additional  parameters  for  phase  rotation.  These  parameters  are 
determined  at  the  receiver  in  the  sense  of  maximizing  the  minimum  channel  gain 
or  capacity,  and  they  are  sent  back  to  the  transmitter  along  with  the  beamforming 
vectors  through  the  feedback  channel. 

For  practical  implementation,  the  proposed  interpolator  with  quantized  phase 
rotations  is  combined  with  existing  approaches  for  beamforming  vector 
quantization.  The  performance  of  the  proposed  beamforming  scheme  is  found  to 
outperform  comparable  beamforming  techniques  in  terms  of  bit  error  rate  (BER) 
and  capacity  and  yet  requires  less  channel  state  information. 
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Abstract — Transmit  beamforming  and  receive  combining 
can  be  simply  implemented  over  frequency  selective  channels 
by  using  multiple-input  multiple-output  (MIMO)  technology 
combined  with  orthogonal  frequency  division  multiplexing 
(OFDM)  in  what  is  known  as  MIMO-OFDM.  Generally  this 
requires  knowledge  of  the  transmit  beamforming  vectors  for 
every  tone  at  the  transmitter.  In  systems  without  reciprocity 
between  the  uplink  and  the  downlink  channels,  this  requires 
informing  the  transmitter  about  all  the  beamforming  vectors 
through  a  feedback  control  channel.  To  reduce  the  feedback 
required,  this  paper  proposes  to  exploit  channel  coherence  to 
send  back  only  a  fraction  of  the  beamforming  vectors.  The 
coherence  of  the  transmit  beamforming  vectors  is  established. 
A  new  spherical  interpolator,  with  an  additional  parameter  for 
phase  rotation,  is  proposed  to  fill  in  the  gaps  between  known 
vectors.  The  parameters  of  the  interpolator  are  optimized  in 
the  sense  of  maximizing  mutual  information  and  are  sent 
back  to  the  transmitter  along  with  the  beamforming  vectors 
through  the  feedback  channel.  Simulation  results  show  that 
the  proposed  scheme  improves  performance  and  reduces 
feedback  requirements  compared  to  existing  beamforming 
techniques. 

i.  Introduction 

The  spatial  dimension  of  multiple-input  multiple-output 
(MIMO)  systems  can  be  used  to  mitigate  signal-level 
fluctuations  in  fading  channels  through  thoughtful  exploita¬ 
tion  of  antenna  diversity.  In  a  MIMO  system,  diversity 
can  be  obtained  by  employing  space-time  codes  [1]— [3] 
without  transmit  channel  knowledge  or  alternatively  by 
using  the  channel  state  information  (CSI)  at  the  transmitter 
using  transmit  beamforming  and  receive  combining  [4]— [8] . 
Compared  with  space-time  coding,  beamforming  and  com¬ 
bining  achieve  full  diversity  order,  they  work  with  any 
number  of  transmit  and  receive  antennas,  and  they  achieve 
additional  array  gain. 

Diversity  techniques  proposed  for  narrowband  channels 
can  be  easily  extended  to  frequency  selective  MIMO  chan¬ 
nels  by  employing  orthogonal  frequency  division  multiplex¬ 
ing  (OFDM).  In  a  MIMO-OFDM  system  with  channel  state 
information  [9],  transmit  beamforming  with  receive  com¬ 
bining  is  performed  independently  for  each  subcarrier  [9], 
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[10].  This  approach,  however,  requires  the  knowledge  about 
the  transmit  beamforming  vectors  at  the  transmitter.  In 
non-reciprocal  channels,  this  necessitates  that  the  receiver 
informs  the  transmitter  about  the  beamforming  vectors 
for  all  subcarriers  through  a  feedback  control  channel.  A 
practical  solution  for  reducing  the  amount  of  feedback 
per  beamforming  vector  is  to  use  a  codebook  designed 
for  quantization  of  the  transmit  beamforming  vectors  (see 
e.g.,  [1 1]— [13]).  The  feedback  requirements,  though,  still 
increase  in  proportion  to  the  number  of  subcarriers. 

In  an  attempt  to  reduce  the  feedback  requirements, 
in  this  paper  we  propose  a  new  approach  to  transmit 
beamforming  that  combines  limited  feedback  of  beamform¬ 
ing  information  and  interpolation  of  beamforming  vectors. 
In  OFDM,  the  number  of  subcarriers  is  designed  to  be 
much  larger  than  the  length  of  the  discrete  time  domain 
channel  impulse  response  to  increase  spectral  efficiency 
(see  e.g.,  [14],  [15]).  This  causes  significant  correlation 
between  neighboring  OFDM  subchannels  which  results  in 
high  correlation  between  beamforming  vectors.  We  analyze 
the  correlation  between  beamforming  vectors  and  use  this 
result  to  reduce  the  amount  of  feedback  information.  We 
propose  to  subsample  the  space  of  beamforming  vectors 
and  send  back  only  the  beamforming  vectors  for  select 
subcarriers  to  the  transmitter.  The  transmitter  in  turn  de¬ 
termines  the  missing  beamforming  vectors  through  a  smart 
interpolation  of  the  conveyed  beamforming  vectors. 

In  this  paper,  we  assume  that  the  transmit  power  is 
identically  assigned  to  all  subcarriers.  The  beamforming 
vectors  have  unit  norm  by  the  transmit  power  constraint, 
thus  the  outputs  of  the  interpolation  problem  should  main¬ 
tain  the  unit  norm  constraint.  This  is  known  as  spherical 
interpolation  [16]— [18].  Unfortunately,  it  is  not  easy  to 
directly  use  spherical  interpolation  algorithms  to  interpolate 
beamforming  vectors  because  the  optimal  beamforming 
vectors  are  not  unique.  As  an  alternative,  we  propose  a 
new  spherical  linear  interpolator  with  additional  parameters 
for  phase  rotation,  which  consider  the  nonuniqueness  of 
the  optimal  beamforming  vector.  Our  interpolator  can  be 
used  to  reduce  the  feedback  requirements  for  transmit 
beamforming  and  can  be  used  to  reduce  complexity  at  the 
receiver  by  permitting  the  singular  value  decomposition  to 
be  calculated  for  only  a  fraction  of  tones.  To  maximize 


Fig.  1.  Block  diagram  of  a  MIMO-OFDM  system  with  Mt  transmit 
antennas,  Mr  receive  antennas,  and  N  subcarriers. 


the  performance  of  the  proposed  interpolator,  the  phase 
parameters  are  determined  in  the  sense  of  maximizing 
the  capacity.  In  fact,  it  is  not  trivial  to  find  the  optimal 
solution  from  the  cost  function  maximizing  the  capac¬ 
ity.  Fortunately,  we  can  use  a  numerical  grid  search  or 
obtain  a  closed  form  solution  through  an  approximation 
of  the  cost  function.  Computer  simulations  show  that  the 
proposed  scheme  outperforms  existing  diversity  techniques 
with  comparable  feedback  information  and  performs  close 
to  the  ideal  beamforming  with  feedback  of  all  beamforming 
vectors. 

For  practical  implementation,  the  proposed  interpolation 
method  can  be  combined  with  beamformer  quantization 
techniques  illustrated  in  [1 1]— [13].  In  this  case,  the  quan¬ 
tized  beamforming  vectors  are  used  to  find  the  optimal 
phases  and  perform  the  interpolation.  We  defer  the  details 
of  quantizing  the  beamforming  vector  and  the  phase  infor¬ 
mation  are  discussed  in  [19]. 

This  paper  is  organized  as  follows.  Section  2  intro¬ 
duces  the  MIMO-OFDM  system  with  beamforming  and 
combining.  We  presents  the  proposed  interpolation  based 
beamforming  scheme  and  the  phase  optimization  for  the 
proposed  spherical  interpolator  in  Section  3.  The  proposed 
beamforming  method  is  compared  with  existing  diversity 
techniques  in  terms  of  bit  error  rate  (BER)  and  capacity  in 
Section  4.  Finally,  we  provide  some  conclusions  in  Section 
5. 


the  subcarrier  fc,  the  received  signal  after  processing  with 
the  combining  vector  z(k)  =  [zi(k),  22 (fc),  •  •  •  ,  ZMr{k)]T 
can  be  expressed  as 

r(k)  =  zH  (k){H(k)w(k)s(k)  +  n(fc)},  1  <  k  <  N  (1) 


where  n (k)  is  the  Mr -dimensional  noise  vector  whose 
entries  have  the  independent  and  identically  distributed 
(i.i.d.)  complex  Gaussian  distribution  with  zero  mean  and 
variance  No.  We  assume  that  the  power  is  allocated  equally 
across  all  subcarriers,  thus  E[\s(k)\2]  =  Ss  is  a  constant 
and  ||w(fc)||  =  1  (where  ||(-)||  means  2-norm  of  (•))  to 
maintain  the  overall  power  constraint. 

Since  the  signal  model  in  (1)  is  identical  to  that  of  a 
narrowband  MIMO  system,  w (k)  and  z(k)  can  be  chosen 
to  maximize  the  signal  to  noise  ratio  (SNR)  using  the 
results  for  narrowband  MIMO  systems  [4],  [6].  Without 
loss  of  generality,  we  can  fix  ||z(fc)||  =  1.  Then  the  SNR 
for  subcarrier  k  can  be  written  as 


ftWBWlf.4  (2) 

iVo  iVo 

where  T(k)  =  \zH  (k)H(k)w(k)\2  is  the  effective  channel 
gain. 

Given  w(fc),  it  is  possible  to  show  that  the  SNR  max¬ 
imizing  solution  uses  maximum  ratio  combining  (MRC) 
with 

H(Jfe)w(fc) 


Z  (k)  = 


l|H(*)w(fc)ir 


(3) 


On  the  other  hand,  given  z(k),  the  effective  channel  gain 
T (k)  can  be  maximized  through  maximum  ratio  transmis¬ 
sion  (MRT).  With  MRT  and  MRC,  w(fc)  is  simply  the  dom¬ 
inant  right  singular  vector  of  H(fc)  corresponding  to  the 
largest  singular  vector  of  H (fc)  [6].  The  MRT/MRC  scheme 
obtains  full  diversity  order  in  Rayleigh  fading  channels  and 
achieves  the  full  array  gain  available.  Unfortunately,  this 
improved  performance  comes  at  the  expense  of  knowledge 
of  the  beamforming  vectors  {w (fc),  k  =  1,  2,  •  •  •  ,  TV}.  The 
extensive  feedback  required  for  the  optimal  MRT/MRC 
solution  motivates  developing  new  methods  that  achieve 
near  maximum  SNR  performance  but  with  more  realistic 
feedback  requirements. 


2.  System  Overview 

Fig.  1  illustrates  a  MIMO-OFDM  system  with  trans¬ 
mit  beamforming  and  receive  combining  with  Mt  trans¬ 
mit  antennas,  Mr  receive  antennas,  and  N  subcarriers. 
During  one  OFDM  symbol  period,  the  scalar  symbol 
s(k)  for  the  k- th  subcarrier  is  mapped  into  the  trans¬ 
mit  antennas  using  the  beamforming  vector  w (fc)  = 
[wi(k),  W2(k),  •  •  •  ,WMt{k)]T.  Assuming  the  sampled  im¬ 
pulse  response  of  the  channel  is  shorter  than  the  cyclic 
prefix  (CP),  the  channel  for  the  k- th  subcarrier  can  be  de¬ 
scribed  by  a  Mr-by-Mt  channel  matrix  H(fc)  whose  entries 
represent  the  channel  gains  experienced  by  subcarrier  k.  At 


3.  INTERPOFATION  BASED  BEAMFORMING 
FOR  MIMO-OFDM 

In  this  section,  we  analyze  the  correlation  between  sub¬ 
channels  of  MIMO-OFDM  and  show  that  the  correlation 
between  the  beamforming  vectors  is  similar  to  the  subchan¬ 
nel  correlation  in  a  multiple- input  single-output  (MISO) 
system.  Next  we  use  these  results  to  motivate  subsampling 
the  beamforming  vectors  followed  by  beamformer  interpo¬ 
lation.  We  introduce  a  beamformer  interpolator  based  on 
spherical  linear  interpolator  but  with  additional  parameters 
for  phase  rotation,  which  takes  the  invariance  of  the  optimal 
beamforming  vectors  into  account.  Finally,  we  optimize 


the  phase  rotation  parameters  in  the  sense  of  maximizing 
mutual  information. 


3.1.  Channel  Correlation  Between  Subcarriers 

Let  us  define  the  L-dimensional  vector  gij  = 
[^j(O), ^j(2),  •  •  •  1)]T  which  denotes  the  time- 

domain  channel  impulse  response  with  length  L  between 
the  transmit  antenna  j  and  the  receiver  antenna  i.  Suppose 
that  {gi,j(n)}  are  i.i.d.  complex  Gaussian  with  zero  mean 
and  unit  variance.  For  notational  convenience,  we  define 
the  (MtMr)~ dimensional  vector  representing  the  channel 
gains  for  subcarrier  k  as 

h(Jfe)  =  vec(H(Jfe))  (4) 

where  vec(X)  means  the  vector  obtained  by  stacking  the 
columns  of  X  on  top  of  each  other  and  1  <  k  <  N. 
The  channel  vector  h(fc)  is  obtained  by  the  discrete  Fourier 
transform  (DFT)  of  time-domain  channel  impulse  responses 
and  hence  denoted  by 

h(fc)  =  Grffe  (5) 


where  G  =  [gM,  •  •  •  ,gMr,i»Si,2,  •  •  •  ,gMr>2,  •  •  •  ,g Mr,Mt] 

is  the  L-by-MtMr  matrix  composed  of  all  time-domain 
channel  impulse  responses  between  transmit  and  receive 
antenna  pairs,  and  is  the  first  L  elements  of  the  k- th 
column  of  the  DFT  matrix  which  is  given  by 


ffe  =  X[l,e~j2nk/N ,e~j2n(-2k)/N ,  ■  ■  ■  ,e-i2T'(L~1)k/N]T , 

N  (6) 
To  measure  the  correlation  between  the  subchannel  k  and 
the  subchannel  (fc  +  d),  a  new  real  variable  ph(d)  is  defined 
as 

E[\hH  (k  +  d)h(k)\2} 


Ph(d)  = 


mmw4 


(7) 


where  0  <  Ph(d)  <  1.  Suppose  that  the  channel 

gains  between  different  antenna  pairs  are  independent,  i.e. 
E[GHG]  =  LIl.  Then  ph(d)  is  expressed  as 


Ph(d) 


t»+dP(k)tk+d 
E[\hH(k)h(k)\2] 
iy2ffp(o)fd 
L2MtMr{MtMr  + 1) 


(8) 


where  P(fc)  = 

mathematical 
as  the  L  x 

r(MtMr+L)MtMr 
L  N  ’ 


=  £[G*GTffcf^G*GT].  By  some 
manipulations,  P(0)  is  denoted 
L  circulant  matrix  which  has 

(MtMr)2  (MtMr)2 

N  ’  '  '  '  ’  N 


as 


its  first 


row. 

In  a  similar  manner  to  the  subchannel  correlation,  the 
correlation  between  beamforming  vectors  is  defined  as 


pw(d)  =  E[\wH(k  +  d)w(fc)|2]  (9) 


Fig.  2.  Comparison  between  p\1(d)  with  one  receive  antenna  and  pw(d) 
with  two  receive  antennas  when  Mt  =  4,  N  =  64,  and  L  =  6  or  12. 


H (fc),  it  is  not  easy  to  express  w (k)  in  terms  of  h(fc)  in  (4). 
This  makes  it  difficult  to  obtain  a  closed-form  expression 
for  pw(d). 

As  an  alternative,  the  correlation  between  beamforming 
vectors  are  examined  through  computer  simulations.  Fig.  2 
shows  pw(d)  when  Mt  =  4,  Mr  =  2,  TV  =  64,  and 
L  =  6  or  12.  For  comparison,  ph(d)  given  by  (8)  were 
plotted  for  Mt  =  4,  Mr  =  1,  N  =  64,  and  L  =  6  or 
12.  Of  particular  interest  is  that  the  correlation  between 
beamforming  vectors  pw(d)  is  very  similar  to  the  channel 
correlation  ph(d)  with  one  receive  antenna.  In  [12],  it  was 
shown  that  the  distribution  of  the  optimal  beamforming 
vector  is  independent  of  the  number  of  receive  antennas. 
Thus,  we  can  use  the  coherence  bandwidth  predicted  by 
Ph(d)  with  one  receive  antenna  to  estimate  the  coherence 
of  the  beamforming  vectors. 

3.2.  Proposed  Interpolation  Based  Beamforming  Method 

A  simple  beamforming  method  for  using  the  correlation 
of  the  beamforming  vectors  is  to  combine  the  neighboring 
subcarriers  into  a  cluster  and  use  the  beamforming  vector 
corresponding  to  the  center  subcarrier  in  the  cluster.  This 
method  will  be  referred  to  as  clustering  1 .  If  we  combines 
K  subcarriers  into  one  cluster,  the  amount  of  feedback 
information  is  reduced  to  1/K .  However,  as  the  cluster  size 
K  increases,  the  beamforming  performance  is  significantly 
degraded  because  of  the  distortion  experienced  by  the 
subcarriers  near  the  cluster  boundary. 

To  mitigate  the  performance  degradation  by  clustering, 
we  consider  a  new  beamforming  method  that  reconstructs 
the  beamforming  vectors  for  all  subcarriers  through  interpo¬ 
lation.  In  the  proposed  scheme,  the  receiver  subsamples  the 


where  ||w(fc)||  1  for  all  k.  Since  the  optimal  beamform  ijn  [20],  a  clustering  scheme  based  on  the  channel  correlation  was 

ing  vector  w(fc)  is  the  dominant  right  singular  vector  of  proposed  to  reduce  the  feedback  information  for  adaptive  modulation. 


subcarriers,  i.e.  selects  a  fraction  of  OFDM  subcarriers,  and 
evaluates  the  optimal  beamforming  vectors  for  the  selected 
subcarriers.  Suppose  that  the  channel  matrices  {H(fc),  1  < 
k  <  N}  are  known  to  the  receiver  and  that  N  is  divided 
by  the  subsampling  rate  K.  The  receiver  evaluates  the 
optimal  beamforming  vectors  for  the  selected  subcarriers 
{w(l),  w (K  +  1),  w(2 K  +  1), • • •  ,  w (N  —  K  +  1)}  and 
sends  them  back  to  the  transmitter.  Then  the  transmitter 
evaluates  the  beamforming  vectors  for  all  the  subcarriers  by 
interpolating  the  conveyed  vectors.  Since  the  beamforming 
vectors  have  unit  norm  for  power  constraint  reasons,  it  is 
natural  to  employ  the  spherical  interpolation  algorithms 
in  [16]— [18].  Unfortunately,  these  algorithms  cannot  be 
applied  to  our  approach  for  the  following  reason.  Recall 
the  effective  channel  gain  T(k)  =  ||H(&)w(&)||2  in  (2). 
When  w (k)  is  the  optimal  beamforming  vector  maximiz¬ 
ing  the  effective  channel  gain,  e^w(k)  also  maximizes 
the  effective  channel  gain.  In  other  words,  the  optimal 
beamforming  vector  is  not  a  unique  point  but  a  line  on 
the  unit  sphere.  As  a  consequence,  algorithms  that  are 
used  to  compute  the  optimal  beamforming  vector,  e.g.  the 
singular  value  decomposition,  typically  choose  the  phase 
to  force  the  first  coefficient  of  each  w (k)  to  be  real.  The 
phase,  however,  has  a  dramatic  impact  on  the  resulting 
interpolations.  This  makes  it  difficult  to  directly  apply  the 
spherical  interpolation  algorithms  to  the  interpolation  of 
beamforming  vectors. 

To  overcome  this  difficulty,  we  propose  a  new  algorithm 
through  modification  of  the  interpolator  in  [18].  Given 
{w(l),  w (K  +  1),  w(2 K  +  1),  •  •  •  ,  w (TV  —  K  +  1)},  the 
proposed  interpolator  is  expressed  as 

w(lK  +  k;0t)  = 

(1  -  cfc)w(Z-ff  +  1)  +  cfc{e-?eiw((7  +  1  )K_  +  1)} 

||(1  -  ck)w(lK  +  1)  +  ck{e^w((l  +  1  )K  +  1)}|| 

(10) 

where  C&  =  (k  —  l)/K  is  the  linear  weight  value, 
w (TV  +  1)  =  w(l),  1  <  k  <  K,  Si  is  a  parameter 
for  phase  rotation  with  0  <  /  <  N/K  —  1.  Note  that 
{w(&) ,7V  —  K -\-l  <  k  <  TV}  are  obtained  by  w(7V  —  K  + 
1)  and  e^eN/K~1w(N  +  1)  (equivalently  w(l)). 

While  the  spherical  interpolators  in  [16]— [1 8]  only  utilizes 
w(lK  +  l)  and  w((/  +  1)K  +  1),  the  proposed  interpolator 
evaluates  the  beamforming  vector  from  w (IK  +  1)  and 
ei°l  w((Z  +  1)K  + 1).  Essentially  the  role  of  Oi  is  to  remove 
the  distortion  caused  by  the  arbitrary  phase  rotation  of  the 
optimal  beamforming  vectors. 

3.3.  Phase  Optimization  for  the  Proposed  Interpolator 

In  this  subsection,  the  proposed  interpolator  is  optimized 
by  determining  the  phases  {6^,0  <  l  <  N/K  —  1} 
in  the  sense  of  maximizing  the  capacity.  The  determined 


phase  parameters  are  conveyed  back  along  with  the  selected 
beamforming  vectors  to  the  transmitter. 

In  (10),  Oi  is  only  used  for  computing  {w (IK  +  fc),  1  < 
k  <  K},  and  thus  the  optimal  6i  maximizing  the  capacity 
can  be  found  by 

A,  f1 ,  iiH(^+fc)w(^+M)n 

0i  =  argmax^loga  jl  + - — - > 

(ID 

where  0  <  /  <  N/K  —  l.  Due  to  the  normalization  factor  in 
w (IK+k;  0 ),  it  is  not  easy  to  get  a  closed-form  solution  for 
(11).  Instead,  we  use  a  numerical  grid  search  by  modifying 
(11)  as 


K 


0i  =  argmaxy^log2  l  1 

60  A--1  '• 


||H(ZA:  +  fc)w(^  +  A;;6>)||: 
No 


'}■ 


(12) 

where  ©  =  {0,  •  ,  2(p^1)71'},  an(j  p  ^  number 

of  quantized  levels  which  determines  the  performance  and 
complexity  of  the  search. 

As  an  alternative  to  the  grid  search,  (11)  is  simplified 
by  considering  the  average  channel  capacity.  In  general, 
w (IK  +  K/2  +  1)  experiences  the  largest  distortion  by 
interpolation,  because  the  subcarrier  (IK  +  K/2  + 1)  is  the 
furthest  from  the  subcarriers  (IK  P  1)  and  ((/  +  1)K  +  1) 
with  feedback  of  the  beamforming  vectors.  Using  this 
observation,  (11)  can  be  approximated  to  maximize  the 
capacity  of  the  subcarrier  (IK  +  K/2  +  1).  Then,  9i  is 
approximately  obtained  by 


Oi  =  arg  max  log2  <  1  + 

0  l 

||H (IK  +  K/2  +  l)w (IK  +  K/2  +  1;  0)\\2 

N0 

By  the  concavity  of  log  function,  this  equation  is  equiva¬ 
lently  rewritten  as 


Oi  =  arg  max  ||H(/iV  +  K/2  +  l)w (IK  +  K/2  +  1;  ( 
6 


=  arg  max  - 
6 


(wi 


pj6\ 


W2)^R(wi 


pjO 


W2) 


||  wi  +  e-7'6>w2||2 


(14) 


where  wi  =  w (IK  +  1),  w2  =  w((Z  +  1  )K  + 1),  and  R  = 
H h(IK  +  K/2  +  1)H (IK  +  K/2  +  1).  By  differentiating 
the  cost  function  in  (14)  with  respect  to  0 ,  we  get 


Imag^e7*9  +  a2)  =  0  (15) 


where  Imag(-)  means  the  imaginary  part  of  (•),  ai  = 
(wfRw2)(wfwi  +  wfw2)  —  (w{iRwi)(wfw2)  — 
(w2ffRw2)(wfw2),  and  a2  =  2(wfRw2)(wfwi).  It 
can  be  shown  that  |aq|  >  |Imag(a2)|  so  that  there 
always  exists  6  satisfying  (15).  Suppose  that  sin-1  (a?)  E 
[— 0.57T,  0.57r)  and  let  denote  oti  =  |a^|eJ<^.  Then  the 


Fig.  3.  Capacity  without  quantization  of  beamforming  vectors  when 
Mt  =  4,  Mr  =  2,  N  =  64,  and  K  =  8.  Grid  search  in  (12)  and 
closed-form  solution  in  (14)  were  used. 


Fig.  4.  Comparison  of  BER  performance  when  Mt  =  4,  Mr  =  2, 
N  =  64,  K  =  8,  and  channel  coding  was  used. 


solutions  of  (15)  are  given  by 


#2 


—01  +  e 

f  —0i  +  7T  —  e  if  e  >  0 
\  —01  —7 r  —  e  if  e  <  0 


(16) 

(17) 


where  e  =  sin  1  1^|  sin  02^  •  Also,  the  derivative  of 

(15)  is  denoted  as 

Mjj{a l^9  +  —  c^2 )  —  OL\e~$0}  =  —2  Real(aieJ'0) 

(18) 

where  Real(-)  means  the  real  part  of  (•).  Since  (18)  gives 
the  sign  of  the  second  derivative  of  the  cost  function, 
the  optimal  6i  maximizing  the  cost  function  is  one  of 
the  solutions  $i  and  $2  such  that  —2  Real(a'ieJ'0*)  is 
negative.  This  solution  will  be  referred  to  as  the  closed- 
form  solution. 


4.  Simulation  Results 

Through  computer  simulations,  we  compare  the  pro¬ 
posed  beamforming  method  with  antenna  selection  di¬ 
versity  [7],  the  simple  clustering  method,  and  the  ideal 
beamforming  with  feedback  of  all  beamforming  vectors.  In 
the  simulation,  we  consider  a  MIMO-OFDM  system  with 
parameters:  Mt  =  4,  Mr  =  2,  N  =  64,  K  =  8.  We 
assumed  that  each  discrete-time  channel  impulse  response 
had  6  taps  (. L  =  6)  with  a  uniform  power  delay  profile 
and  i.i.d.  complex  Gaussian  distribution  with  zero  mean 
and  variance  1  / L\  the  channels  between  different  transmit 
and  receiver  antenna  pairs  were  independent;  n (k)  was 
i.i.d.  complex  Gaussian  with  zero  mean;  the  feedback 
channel  had  no  delay  and  no  transmission  error;  receiver 
used  MRC  with  perfect  channel  knowledge;  quadrature 
phase  shift  keying  (QPSK)  modulation  was  used  for  BER 


simulations;  the  transmit  power  was  identically  assigned  to 
all  the  OFDM  subcarriers.  Every  point  of  the  simulation 
results  was  obtained  by  averaging  over  more  than  5000 
independent  realizations  of  channel  and  noise. 

Fig.  3  presents  the  channel  capacity  obtained  by  various 
beamforming  techniques.  In  the  figure,  ’’grid  search”  and 
’’closed-form  solution”  describe  the  proposed  interpolators 
using  the  grid  search  in  (12)  and  the  closed- form  solution 
in  (14),  respectively.  The  proposed  grid  search  and  closed- 
form  solution  had  higher  capacity  than  the  clustering  and 
antenna  selection,  but  showed  slight  loss  compared  to 
the  ideal  beamforming.  The  proposed  methods  performed 
within  0.3  dB  of  the  ideal  beamforming,  while  the  cluster¬ 
ing  had  about  1  dB  loss  at  capacity=4  bps/Hz. 

To  further  evaluate  the  performance  of  the  proposed 
scheme,  the  BER  values  with  channel  coding  and  in¬ 
terleaving  are  presented  in  Fig.  4.  For  channel  coding, 
we  used  a  convolutional  code  with  generator  polynomials 
go  =  1338  and  g\  =  171s  with  coding  rate  1/2,  along 
with  the  interleaver  and  deinterleaver  defined  in  [15],  and 
soft  Viterbi  decoding.  The  frame  length  was  30  OFDM 
symbols  and  each  OFDM  symbol  transmitted  64  bits.  The 
channel  was  assumed  to  be  fixed  for  a  frame  and  randomly 
varied  between  frames.  As  expected  in  the  capacity  result, 
the  proposed  closed-form  solution  performed  comparable 
to  the  grid  search  with  P  =  4,  and  exhibited  0.8  dB  gain 
over  the  clustering  and  2.2  dB  gain  over  selection  diversity 
at  BER=10-4.  As  compared  with  ideal  beamforming,  the 
Eb /No  loss  of  the  proposed  was  within  0.2  dB. 

5.  Conclusions 

A  new  transmit  beamforming  scheme  for  MIMO-OFDM 
was  proposed.  Using  the  correlation  between  neighboring 
beamforming  vectors,  the  proposed  scheme  reduces  the 
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feedback  requirements  for  a  closed-loop  MIMO-OFDM 
system.  A  new  spherical  interpolator  was  proposed  to  inter¬ 
polate  beamforming  vectors  and  the  proposed  interpolator 
can  be  optimized  by  adjusting  the  parameters  for  phase 
rotation  in  terms  of  an  effective  channel  gain  or  capacity. 
Through  computer  simulations,  it  was  shown  that  the  pro¬ 
posed  method  has  considerable  benefits  in  performance  and 
the  amount  of  feedback  information  compared  to  existing 
diversity  techniques.  Our  work  proposed  herein  is  limited  to 
a  MIMO  system  that  transmits  a  single  data  stream  per  sub¬ 
carrier.  Extending  this  approach  to  the  case  where  multiple 
beamforming  vectors,  combining  with  spatial  multiplexing 
as  for  example  in  [22],  remains  a  future  research  topic. 
Accounting  for  imprecision  of  channel  state  information, 
as  in  [10],  is  also  of  interest. 
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Transmit  Diversity  in  MIMO  Systems 


■  No  channel  state  information  (CSI)  at  Transmitter 

■  Space-time  codes  [Tarokh  et  al.  1998] 

■  Full  CSI  at  Transmitter 

■  Antenna  selection  [Wittneben  et  al.  1 994],  [Heath  et  al.  1998] 

■  Maximum  ratio  transmission  [Lo  1999],  [Andersen  2000] 

■  Partial  CSI 

■  Adaptive  beamforming  [Xia  et.  al.  2004] 
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Beamforming  and  Combining 


( with  limited  data  rate )  Complex 

_ number 

■  Covert  MIMO  to  SISO  V  =  z*Hws  +  z*n 

■  T ransmit  beamforming  vector  w 

■  Receive  combining  vector  z 
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Max  SNR  Solution  (MRT/MRC) 


SNR  is  proportional  to  ||z*Hw 


Max  SNR  solution  is 


w  =  max 


w 


=1 


Hw 


z  =  Hw 

Solution  is  not  unique:  if  w  is  optimal,  so  is  weJ0 


■  Other  solutions  possible  (constrain  w) 

■  Selection  diversity 

■  Equal  gain  combining  /  transmission 
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Beamforming  in  MIMO-OFDM 


s(l) 

s(N) 


r(l) 

r(N) 
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Signal  Model  for  One  Subcarrier 


■  Equivalent  model  for  subcarrier  k 

(identical  with  a  narrowband  MIMO  case) 


4*) 


„  k  ,  (k) 


V'' 


:v 


h 


Mr ,  Mt 


(k) 


nM,  ( k ) 


Noise 


nx(k)  Zl{k) 


r{k) 


'M, 


(k) 


r(k )  =  zH  (k){H.(k)w(k)s(k)  +  n(fc)} 
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Problem  Summary 


J\T _  -1 

■  Beamforming  in  MIMO-OFDM  requires  {w(n)}n=0 

■  Feedback  requirements  <x  Number  of  subcarriers 

■  How  can  we  reduce  the  number  of  vectors  fed  back? 

■  Exploit  correlation  between  vectors 

■  Send  back  fraction  of  vectors 

■  Use  “smart”  interpolation 

■  How  can  limit  the  feedback  for  each  vector? 

■  Use  quantized  beamforming  [Love  et.  al.  2003] 
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Correlation  Between  Subcarriers 


♦ 

♦ 


ph(d)\  Mf=4,  /W=1,  W=64,  K=  8,  L={ 6,12} 

Pw(cO:  Mp4,  M=2,  N=64,  K=8,  L={6,12} 

subsampling  rate  number  of  channel  taps 
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Clustering  of  Subcarriers 


■  Clustering  (K=  5) 

cluster  1  cluster  2  cluster  3  cluster  4 


(  ^ 


w(3)  w(8)  w(13)  w(18)  subcarriers 


feedback 


■  Disadvantages 

■  Performance  degradation  in  cluster  boundary 

■  Cluster  size  (or  feedback  reduction)  is  limited 
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Proposed  Beamforming  Method 


Subsampling  of  beamforming  vectors 


4 

/ 


k 

\ 

\ 

\ 


w(l) 

feedback  | 


i 

i 

w(K+l)  w(2K  +  1) 


subcarriers 


w (N  -  K  +  1) 


1 


1 


i 


□  Reconstruction  of  beamforming  vectors 

♦  Transmit  power  constraints  force  ||w(fc)||  =  1 
- *  Spherical  interpolation 
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Conventional  Spherical  Interpolation 


□  Spherical  averaging  [Watson  1983],  [Buss  et  al.  2001] 

IIELoM 

P 

+  p  is  the  order,  ||v||  =  1,  bj  =  1,  and  bj  >  0 

i= 0 

□  Nonuniqueness  of  optimal  beamforming  vectors 

w(fc)  is  optimal  * — ►  e^ew(k)  is  optimal 

♦  Random  phase  rotation  0  has  significant  impact  on  interpolations 

♦  Performance  degradation  of  spherical  interpolators 
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Proposed  Interpolator 


■  Spherical  Linear  Interpolator  with  phase  rotation 


w(fe,6»i) 


(1  -Cfc)w(l)  +cfc{ 

eJ'^i 

*r(K+  1)} 

11(1  -Cfc)w(l)  +Cfc{ 

«r(K+l)}|| 

1  <  k  <  K 


and  0 1 


is  a  parameter  for  phase  rotation. 


□  Optimization  of  phase  rotation  parameters 

♦  Maximize  the  diversity  gain 

♦  Maximize  the  mutual  information 
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Optimization  of  Phase  Rotation 


■  Maximizing  the  minimum  channel  gain  (diversity) 

9\  =  arg  ^  max  min{||H(fc)w(fc,  0)||2,  1  <  k  <  K} 

■  Difficult  to  get  a  closed-form  solution 


■  Grid  search 


01  =  arg  max min{||H(fc)w(fc, 0)|| 


1  <  k  <  K} 


■  ©  =  {o, 

levels 


27 T  47T 
~P’~P' 


2  (P-  1)7 r 
P 


and  P  is  the  number  of  quantized 
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Alternative  Solution 


■  Observation 

■  Subcarrier  (K/2+1 )  suffers  from  the  largest  distortion 

■  Subcarrier  (K/2+1 )  has  the  worst  average  channel  gain 

■  Approximation  of  cost  function 

&i  =  arg  max  ||H(iC/2  + l)w(K/2  + 1,6>)||2 

\\H(K/2  +  l){w(l)  +  e^w  (K  +  1)}|| 

=  arg  max  - - - - — — ^ — - - - — 

6e[  0,27t)  w(1)  +  e^w(K  +  1) 
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Closed-Form  Solution 


■  Differentiation  with  respect  to  0 

Imag^ie^  +  a 2)  =  0 

■  a.,  and  a2  are  complex  constants  satisfying  iai|  >  |imag(a2)l 

■  Always,  there  exist  two  solutions. 

■  Second  derivative 

— Imag(aie:?61  +  =  — 2Real(aie:?61) 

dO 

■  The  solution  making  the  second  derivative  be  negative  is  the 

optimal  solution.  - *  closed-form  solution 
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Optimization  in  Terms  of  Capacity 


Maximizing  the  capacity 

K  ( 

01  =  arg  max  V  log2  1  + 

6>G[  0,271-)^  V 


|H(fc)w(M)l| 

N0 


Grid  search 


K  ( 

Oi  =  arg  max  £,  log2  1  + 

6e°  k=  1  V 


©  =  {0,^,..,^^} 


||H(fc)w(M)l| 

N0 


P  ’  P  ’ 


P 


■  Closed-form  solution  is  the  same  as  in  diversity  case 
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Simulation  Environments 


■  MIMO  channels 

■  Each  channel  impulse  response  has  6  taps  (L= 6)  with  uniform 
profile  and  follows  i.i.d.  CN(0,ML). 

■  The  channels  between  different  transmit  and  receiver  antennas 
pairs  are  independent. 

■  Noise:  i.i.d.  with  CN(0,A/0) 

■  Simulation  parameters 

■  M,= 4  (#  of  Tx  antennas),  M=  2  (#  of  Rx  antennas) 

■  N= 64  (#  of  subcarriers),  K=  8  (subsampling  rate),  L= 6 

■  Basic  assumptions 

■  MRC  at  the  receiver 

■  No  delay  and  no  transmission  error  in  the  feedback  channel 

■  QPSK,  no  water-filling 
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Simulation  Results  -  Channel  Gain 


■  OFDM  subcarrier  vs.  effective  channel  gain 
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Simulation  Results  -  Uncoded  BER 


■  No  channel  coding 


Eb/No(dB) 
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Simulation  Results  -  Capacity 


■  No  power  control  for  diversity  techniques 
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Simulation  with  Quantization 


■  Quantization  of  beamforming  vectors 

■  By  using  the  codebook  in  [Love  &  Heath  2003] 

■  6  bits  of  feedback  for  each  beamforming  vector 

■  Phase  optimization 

■  Grid  search  with  P= 4  (2  bits  of  feedback  per  phase) 

■  Feedback  requirements 

■  Ideal:  6/V  =  384  bits  /  OFDM  symbol 

■  Clustering:  6N/K  =  48  bits  /  OFDM  symbol 

■  Selection:  2N  =  128  bits  /  OFDM  symbol 

■  Proposed:  8N/K  =  64  bits  /  OFDM  symbol 

■  Orthogonal  STBC  (OSTBC):  none 
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Simulation  Results  -  Channel  Gain 


■  Mt= 4,  Mr 2,  A/=64,  K=  8,  L= 6 
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Simulation  Results  -  BER 


■  Selection 

■  Full  diversity  order 

■  Array  gain  loss 

■  Proposed 

■  Diversity  order 
loss 

■  Additional  array 
gain 
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Simulation  Results  -  BER 


■  BER  degradation  by  channel  prediction  error 


♦  Prediction  error  model 

H(fc)  =  7HO)  +  yi  -72U(fc) 
where 

U(/c)  is  a  prediction  error 


10  1 - 1 - 1 - 1 - 

-20  -15  -10  -5  0 

MSE  of  channel  predictor  (dB) 
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Simulation  Results  -  Capacity 


■  No  power  control 
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Simulation  Results  -  Coded  BER 


1/2  convolutional  code 
Interleaving/deinterleaving 
Frame  length 
■  30  OFDM  symbols 

64  bits  /  OFDM  symbol 
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Conclusions 


■  Summary 

■  Interpolation  for  transmit  beamforming  in  MIMO- 
OFDM  proposed 

■  Relates  to  problem  of  spherical  quantization 

■  Future  work 

■  Proposed  system  does  not  obtain  MIMO  capacity 

■  Develop  “interpolated”  waterfilling  solutions 
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